Temperature Dependence of the Conductivity Sum Rule in the Normal State due to 

Inelastic Scattering 
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We examine the temperature dependence of the optical sum rule in the normal state due to inter- 
actions. To be concrete we adopt a weak coupling approach which uses an electron-boson exchange 
model to describe inelastic scattering of the electrons with a boson, in the Migdal approximation. 
While a number of recent works attribute the temperature dependence in the normal state to that 
which arises in a Sommerfeld expansion, we show that in a wide parameter regime this contribution 
can be quite small. Instead, most of the temperature dependence arises from the zeroth order term 
in the 'expansion', through the temperature dependence of the spectral function, and the inter- 
action parameters contained therein. For low boson frequencies this circumstance causes a linear 
T-dependence in the sum rule. We develop some analytical expressions and understanding of the 
temperature dependence. 

PACS numbers: 



I. INTRODUCTION 

In the last few years much attention has been focussed 
on the anomalous temperature dependence of the opti- 
cal sum ruleiiSiSi^iSiSiLS. In particular, in a number of 
superconductors in the underdoped and optimally doped 
regime the sum rule exhibits a temperature dependence 
which is anomalous with respect to that expected from 
a conventional BCS condensation. The origin of this be- 
haviour has been the subject of a number of theoretical 
explanations2ii2iiiii2iiii4. It is probably fair to say that 
there is no consensus concerning the origin of the ob- 
served anomaly below the superconducting critical tem- 
perature, Tc- 

Perhaps equally intriguing is the normal state tempera- 
ture dependence>i^*iSiiZ-i'''2i Most often the value of the 
sum rule is plotted versus T^, where T is the absolute 
temperature, and a linear relationship is observedi ^i^i^ 
This appears natural as a straightforward Sommer- 
feld expansion produces a finite temperature correction 
which is quadratic in temperature. However, several 
researchersSiiSiiS have noticed that the observed temper- 
ature dependence is too strong compared with what one 
would expect with realistic bandwidths in the cuprate 
materials (for a quantitative analysis of experimental 
data in cuprates see Ref. f2l]) . Moreover, there are clear 
deviations from dependence)^ along with some indi- 
cations of a linear temperature variation of the spectral 
weighljSiSi so that one is left with the question: within 
a Fermi liquid approach, what governs the temperature 
dependence of the optical sum rule in the normal state? 

This is the question we address in this paper. To be 
concrete we use an exchange boson model (Migdal ap- 
proximation for phonons), which is a weak-coupling ap- 
proach. This approach is believed to be applicable to 



phonons, whose possible relevance to the cuprates has 
been questioned in the last few years^^. It is more suspect 
from a theoretical point of view for magnon exchange; 
nonetheless we simply adopt the Migdal approximation, 
with justification provided by some researchers^^. This 
type of calculation has also recently been carried out 
by Karakozov and MaksimoAiS^. We reinforce their (and 
our) main message: a strong temperature dependence of 
the optical sum rule can arise due to interactions. We 
examine in further detail the dependence on the vari- 
ous parameters of the model. In particular, we show 
that the finite width of the band, which is crucial for 
discussing the sum-rule temperature dependence, intro- 
duces new features in the electron self-energy, which are 
usually overlooked in the infinite-bandwidth approxima- 
tion. Moreover, we show that the sum-rule temperature 
dependence cannot be attributed merely to the tempera- 
ture dependence of the imaginary part of the self energy. 
Indeed, the real part enters in a critical way. 

The outline of the paper is as follows. In the next 
section we outline the formalism, both on the real and 
imaginary axes. The former is more appropriate for an- 
alytical work, whereas the latter is ideal for numerical 
work. We show results for a wide variety of parameters, 
which clearly shows how the temperature dependence can 
go from quadratic to linear in temperature, depending on 
the boson frequency scale. We also outline how the Som- 
merfeld contribution can be determined (even in the pres- 
ence of interactions) . It is small always since we use the 
Fermi energy as a dominant energy scale. We also show 
how a large but non-infinite energy scale actually qual- 
itatively changes 'low energy' results. In section III we 
adopt a sequence of approximations to illustrate precisely 
what is required to achieve qualitative and quantitative 
agreement with the numerical results. We find that the 
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real part of the single particle self energy is required for 
a clear understanding of the temperature dependence of 
the sum rule. Finally, we close with a summary. 



II. FORMALISM 



The " single band" sum rulei24iiiiSi2& can be written 



dv Re [a. 
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(1) 



where is the tight-binding dispersion and Uk is the 
single spin momentum distribution function. In the case 
of two dimensions with nearest neighbour hopping only, 
the quantity in braces is given by the negative of the 
expectation value of the kinetic energy 
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We have explored elsewhere that the kinetic energy is 
a good monitor of the sum rule, even in cases beyond 
nearest neighbour hoppingiiSiS^ Therefore in what follows 
we will use Eq. as it is much simpler to calculate for 
a given model than the left-hand-side of Eq. ^ . 

In general, the momentum distribution function is 
given by 
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where G{k,iiVm) is the single electron Matsubara Green 
function with momentum k and frequency iuJm- The 
icUrn are the Fcrmionic Matsubara frequencies, itj„i — 
iTT(2m — 1), where m is an integer, T is the temperature 
(fcs = 1), and /3 = 1/T is the inverse temperature. The 
spectral function, A{k,uj) = —^ImG{k,iuj,n — * uj + iS), is 
related to the analytical continuation of the Matsubara 
Green function. Finally, /(w) = l/(e^'^-|-l) is the Fermi- 
Dirac distribution function. In the absence of interac- 
tions the spectral function is a temperature-independent 
Dirac delta function, so that nt = f{ek — A*) and Eq. (0 
reduces to; 



(4) 



It is then clear that the only temperature dependence of 
the sum rule Q can be due to the "thermal smearing", 
i.e. the effect contained within the Fermi-Dirac distribu- 
tion function. A Sommerfeld expansion then reveals this 
temperature dependence to be quadratic in T/W^^^ 



where W is the electronic bandwidth; this parameter is 
usually small at temperatures of interest (up to 200 ~ 300 
K). However, in the general case of an interacting system 
Eq. 1211 can be written, by means of Eq. Q, as 



K{T) ^2 deduj N{e)eA{e,uj)f{uj) 



(5) 



where we converted the momentum sum to an energy 
integration in Eq. (0). In the presence of interactions, 
the Green function G is given by the Dyson equation. 



G{k,z) = l/[z-ek-^ik,z)], 



(6) 



where z is anywhere in the complex plane, and S(fc, z) 
is the electron self energy. This quantity, and conse- 
quently the spectral function A{e,uj), will be in general 
temperature dependent, and it is this temperature depen- 
dence that will manifest itself in the total kinetic energy 
©(and therefore the conductivity sum) that we wish to 
understandtSZi Indeed, using for simplicity a constant den- 
sity of states in Eq. (01, an application of the Sommerfeld 
expansion produccs^^ 



K(T) = Ko{T) + K 
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where the second term is what we refer to in what fol- 
lows as the Sommerfeld term (with chemical potential 
assumed to be zero, in the middle of the band) and 
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In the non-interacting case H{lo) reduces to 2uo/W, so 
that the T"^ /W dependence of the Sommerfeld term 
is recovered, as stated above. In the interacting 
case instead dH{uj)/dhJ is affected by the quasiparticle 
renormalizationlS. (see Eq. H15|l below) and can also bring 
additional temperature dependence to the factor of 
Eq. 10. The first term Ka{T) in Eq. ^ reduces to the 
constant value K{T = 0) in the absence on interactions, 
see Eq. Q. However, in the general case it is given by 
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Ko{T)^ I dujH{Lo)^ / duj—j deeA(e,uj). 

(9) 

The key observation, also made (independently) by 
Karakozov and Maksimov^'^ is that Eq. Q can have 
a far more significant temperature dependence than the 
Sommerfeld contribution. Moreover, as we shall investi- 
gate here, due to this effect the temperature dependence 
of the kinetic energy can deviate from the non- interacting 
case. 

As mentioned in the introduction, this work will fo- 
cus on "conventional" interactions, modelled after the 
electron-phonon interaction, as described by an Einstein 
oscillator with frequency uje and coupling strength A. In 
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the non-selfconsistent Migdal approximation, the self en- 
ergy is given by: 

= — \{iuJm' - iU!m)Go{k',iuJrn'), (10) 

k' ,m' 



where the non-selfconsistency is reflected in the fact 
that the non-interacting electron Green function is used 
(hence the subscript zero). This can be made self- 
consistent but the changes are minor because of the large 
bandwidth we will always adopt (in other words we are 
not pursuing the possibility of strong temperature depen- 
dence arising from a low electronic energy scale). Here 



A(z) 



(11) 



where z is anywhere in the complex plane. The band- 
width W is included in Eq. H10() since it is customary 
to have a density of states factor in A (which, among 
other things, makes it dimensionless), and so we use the 
'average' density of states, 1/W. One can readily eval- 
uate Eq. (|10() for various tight-binding models and use 
this in Eqs. (|6I3I2|I . However, we simplify further be- 
fore presenting results and advise the reader that such 
details make little difference in our conclusions. To this 
end we assume a constant density of states for energies 
—W/2 < e < W/2. Then the k-summation becomes a 
single integral over energy, e. The result is 
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where N and / are the Bose and Fermi functions, respec- 
tively. Eq. 11211 can be evaluated either at the Matsubara 
frequencies iuim = ii^T{2m — 1), with m an integer, or 
at the real frequencies attained through direct analyti- 
cal continuation, iojm — > w -I- with i5 an infinitesimal 
imaginary part. To obtain reliable numerical results the 
most straightforward approach is on the imaginary axis, 
as we now describe. 



A. Imaginary axis formulation 

First, one evaluates the integrals in Eq. H12|l numer- 
ically. One can show that the result is pure imaginary 
(because of the particle- hole symmetry), and so one can 
define a "renormalization" function through S(iWm) = 
iujmil — Z{iuj„i)), where Z(iujm) is a real-valued function. 
Next one returns to the kinetic energy; the momentum 
dependence (i.e. the e dependence) in Eq. ^ is simple 
enough (see the first line of Eq. Q with e dependence 
coming only from the explicit dependence shown in Eq. 
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FIG. 1: The negative of the kinetic energy vs. temperature 
for a variety of coupling strengths. Note that the magnitude 
of the temperature variation increases considerably with cou- 
pling strength, and, as explained in the text, this cannot be 
attributed to the impact of interactions on the Sommerfeld 
term. 



©) to be done analytically. The result is: 
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This is easily summed provided asymptotics are used 
for the contribution at high frequencies. Eq. H13(l with 
Z{iujrn) obtained from Eq. (|12|1 provides us with our 
numerical results. 

It is difficult to develop an intuition for what temper- 
ature dependence may emerge from various parameter 
choices. In particular we have uncertainty about the 
energy scale of the boson. Before developing some an- 
alytical approximations we therefore show some results. 
Fig. 1 shows the (negative of the) kinetic energy vs. tem- 
perature for several coupfing strengths: A = 0, 0.5, 1 and 
2 as labeled. We have used a consistent set of parame- 
ters for these results: (i) a bandwidth W — 1 eV, and 
(ii) a boson frequency oje = 25 meV. The first most no- 
ticeable trend is that the absolute value of the kinetic 
energy decreases with increasing coupling strength. This 
is expected, as coupling to a boson 'decoheres' the elec- 
trons. While the momentum distribution function retains 
a discontinuity at the Fermi level (see, for example. Fig. 
1 in Ref. 14), the overall distribution is considerably 
more smeared than in the non-interacting case. Thus, 
in the ground state, higher kinetic energy (i.e. lower 
absolute value) states are occupied while low kinetic en- 
ergy states remain unoccupied. This sort of trend with 
coupling strength would be very difficult to observe in 
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FIG. 2: The negative of the kinetic energy vs. temperature 
for a variety of boson frequencies. Note that with decreas- 
ing phonon frequency not only is the magnitude of the tem- 
perature variation increased, but it clearly crosses over from 
quadratic to linear in temperature. 



the experiments, however, as there is no simple way to 
compare two systems whose only difference is the cou- 
pling strength. A more promising property to focus on is 
the temperature dependence, which clearly varies as the 
coupling strength increases. Before doing this we wish 
to make one important remark concerning our parame- 
ter choices. Note that the bandwidth has been chosen to 
be by far the largest energy scale in the problem. This 
is generally the case in 'conventional' metals. Interest- 
ing effects can arise when the bandwidth is not so large 
compared to the other energy scales (such as the boson 
frequency or the temperature itself)?^ '^"^^ However, here 
we wish to demonstrate that even in the 'conventional' 
case a significant temperature dependence of the kinetic 
energy appears due to interactions. 

As we explained above, for the non-interacting case 
the deviation of K{T) from a constant is given by 
the T^/W Sommerfeld contributioniSiai The finite- 
temperature corrections for the interacting cases shown 
in Fig. 1 are clearly not simply quadratic, and further- 
more the magnitude of the variation increases substan- 
tially with coupling strength. Is this due to alterations in 
the Sommerfeld contribution? The answer is most defini- 
tively no. Indeed, even though the Sommerfeld correction 
(to be further described below) changes by at most 65% 
(at 200 K) from A = to A = 2, this gives a correction of 
about 1.7 meV (at 200 K) for A = 2, whereas the devia- 
tion from the zero temperature result is « 20 meV. Thus 
the bulk of the temperature correction arises from other 
sources. 

How does the energy scale of the Einstein frequency lo e 



affect the temperature dependence? In Fig. 2 we show 
the (negative of the) kinetic energy vs. temperature for 
several Einstein frequencies with the same bandwidth as 
before, and with a moderate coupling strength, A = 1. It 
is clear that the energy scale of the boson frequency has 
a profound effect both on the zero temperature value and 
the temperature dependence itself. In particular, for low 
boson frequencies the temperature dependence is linear; 
in addition, as will be made clear below, the temperature 
dependence comes almost entirely from Kq{T), not the 
Sommerfeld term. In contrast, for high boson frequency 
the temperature dependence resembles that of the non- 
interacting case and a significant fraction comes from the 
Sommerfeld term. 



B. Real axis formulation: Sommerfeld contribution 
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FIG. 3: The mass enhancement parameter vs. temperature 
for (a) uje = 2 meV, and (b) uje = 10 meV. Note the (very) 
slight deterioration of the approximate result as the boson 
frequency increases. Also note that for low boson frequency 
the mass enhancement parameter becomes negative for high 
temperatures (but 'high' being defined relative to the boson, 
not the electronic energy scale). Counterintuitively, the result 
for infinite bandwidth (dashed green curve) becomes poorer 
as the boson energy scale decreases. 



To see how this significant temperature variation 
arises, let us go back to the real-axis formulation ©-{Tj) 
of the kinetic energy and let us analyze the Sommerfeld 
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contribution. To obtain the self energy on the real axis 
a straightforward analytical continuation of Eq. I|12|) is 
required, 



'W/2 



de 



W/2 



N{iUE)+f{e) 
Lu + iS + lue — e 
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(14) 



Writing the real and imaginary parts as S = Si + iS2, 
it is easy to see that d'E2{uj)/duj is zero at zero fre- 
quency. Furthermore, = 0) = 0. Defining, as usual, 

9Ei(a))/9cj|i^=o = ~Ao(r), then we obtain. 
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Substitution into the second term of Eq. ^ yields, for 
the Sommerfeld contribution. 
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Here, the zero frequency imaginary part of the self energy 
(denoted simply by E2(T)) is required; for a constant 
density of states the (exact) result is given by the result 
known as being a good approximation for the infinite 
bandwidth case'^^. 
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FIG. 4: The negative of the kinetic energy vs. temperature 
showing the breakdown for the various contributions. The 
solid (red) curves are for small boson frequency, ue = 'i meV, 
with A = 1. The two thin lines show how the zero tempera- 
ture value is modified by the Sommerfeld term by a very small 
amount; the full numerical result is given by the thick red 
curve. A similar contribution is found for the non-interacting 
case (actually, there it is a bit larger) by the dotted (blue) 
curves. However, here the relative efltect is much larger be- 
cause only the Sommerfeld term is responsible for the tem- 
perature variation. Finally, for higher boson frequency, the 
dashed (green) curves illustrate the amount of temperature 
variation due to the Sommerfeld term compared with the rest. 



S2(r) = -7rAcjB/sinh(/3a;_B). 



(17) 



At low temperatures this is zero, while at high temper- 
atures (with respect to the boson frequency, uje), this is 
linear in temperature. 

The real part of the self energy, required to determine 
Ao(T), is more difficult. For loe < W/2 (the physically 
relevant case) we obtain 



Ao(r) = A 47rT^ 
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where /{uje) is the Fermi function. This expression is 
not very enlightening. However note that the first line 
is strongly convergent, so much so that the tan^^ term 
is not really required for most parameter regimes of in- 
terest. The Matsubara sum can then be performed, and 
one obtains for this first term (without the A out front), 
— W£;lm-0'(i -I- UDe), where uje = u}e/{2-kT) and ^'(•z) 
is the diagamma function; this is the result one obtains 



for the infinite bandwidth case. It is positive definite, 
starts at unity at zero temperature (so Ao(T ~ 0) « A) 
and goes to zero as (w^/T)^ at high temperature. The 
other two terms vanish for infinite bandwidth. In fact the 
third term is very small for all temperatures with con- 
ventional parameters, while the second term dominates 
at high temperature. This is clear since it increases (in 
magnitude) linearly with temperature. A plot is shown 
in Fig. 3. Also shown is the approximation 



IIjJe) H ; — . 

' TT W/2 ' 



(19) 



where we have used Eq. (|17|) . This approximation is 
shown by the symbols and is clearly excellent over the en- 
tire temperature range. The important feature apparent 
in Fig. 3 is that the mass enhancement parameter Aq (T) 
changes sign as the temperature increases. This is a 
non-trivial effect; Cappelluti et al.^^ noted a similar phe- 
nomenon when impurity scattering was included (we re- 
mind the reader that at high temperature low frequency 
phonons behave much like static impurities). This be- 
haviour is necessary to include to describe the numerical 
results accurately, as we shall see below. It is truly re- 
markable that at temperatures below 100 K a modest 
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scattering (A = 1) off low frequency phonons [cue = 2 
meV) exhibits a qualitative change (mass enhancement 
changes sign) when an electronic cutoff energy (normally 
taken to be infinite) is introduced which is 2-3 orders of 
magnitude larger (1 eV) than the low energy parameters. 

Once we have exactly determined the quasiparticle 
renormalization parameter Ao(T) we can estimate the 
Sommerfeld contribution. The results are plotted in Fig. 
4, where we compare again the non-interacting case with 
the A = 1 cases at two phonon frequencies. As one can 
see, in the non-interacting case the Sommerfeld term, 
even if small, is the only source of the temperature varia- 
tion of the kinetic energy. When interactions are present 
instead its contribution cannot account for the temper- 
ature variation of K{T) obtained through the numerical 
computation of Eq. (|13|l . Wc then conclude that to un- 
derstand the origin of the strong temperature dependence 
displayed in Figs. 1-2, one must examine Eq. I^. Note 
that a different conclusion was instead drawn in Ref.lH 
where a Hubbard-like interaction between electrons was 
considered. Indeed, in this case the main effect of the 
interaction is to reduce the effective quasiparticle band- 
width 1/Ao, so that as Aq increases the Sommerfeld term 
continues to provide the largest contribution to the tem- 
perature variation of the kinetic energy. 



III. ORIGIN OF THE LINEAR TEMPERATURE 
DEPENDENCE 

A. Imaginary part of the self energy 

The linear temperature dependence, and, indeed, the 
strongest temperature dependence, occurs for low boson 
frequency, so we focus on this parameter regime. To 
put this another way, one should examine temperatures 
which are high compared to the boson frequency, which 
Fig. 2 illustrates is readily attained in the 100 ^ 200 
K range when loe ^ 25 meV. Then it is clear that the 
imaginary part of the self energy at zero frequency in- 
creases linearly with temperature (see Eq. IjlTf) ') and one 
is tempted to associate the linear temperature depen- 
dence observed in K{T) with the linear temperature de- 
pendence in Yj2{T). In fact, Karakozov and Maksimov^^ 
noted that the infinite frequency limit of the self energy 
(as defined by the theory with infinite bandwidth), which 
we will denote by Y,2oo{T), also increases linearly with 
temperature, and they use this quantity. It is given by 



^2oo(T) = -7r^^^coth(^^; 



(20) 



We should be clear here; the physical self energy always 
has an imaginary part which goes to zero beyond the 
electronic bandwidth energy scale. However, the stan- 
dard theor y?^'?^ e liminates this cutoff, and then Eq. H20|l 
provides the imaginary part of the self energy in the in- 
finite frequency limit. 



Either Eq. H20() or (|17() provides the correct high tem- 
perature behaviour for the imaginary part of the self en- 
ergy. However, as we now illustrate, using this part of 
the self energy alone is not sufficient to properly repro- 
duce the numerical result. To examine what the neces- 
sary ingredients are, we proceed in steps, and illustrate 
what approximations to the self energy are sufficient in 
Eq. lO to reproduce accurately the numerical result. A 
first step (following Ref. l25i) is to model the self energy 
with a constant pure imaginary part, which we take to 
be Ti2ao{T). Then the spectral function is given by 



A{t,uj) = 



1 



{T)\ 



(21) 



over all frequency, lo. Insertion into Eq. then allows 
both integrals to be done analytically. If we first assume 
that Eq. lf?T)) holds for all frequencies (w), then Eq. 
results in 



Koa{T) 



where x = 



W 



-AT)\ 



X J— tan 

TT 



1, 



(22) 



Clearly, as |S2oo(T)l decreases com- 
-W/A, as for the non-interacting 
along with the 
5. 



W/2 

pared with W/2, Kq 
case. This result (green, dashed curve) 
numerical result (solid, red curve), is shown in Fig. 
Note that the numerical result is that for Ko{T); this 
differs from the results presented in Figs. 1 - 2 by essen- 
tially the Sommerfeld term, which we leave out of this 
discussion. In fact we have checked that a numerical 
evaluation of Ko{T) based on the exact self energies in- 
serted into Eq. ^ agrees with the imaginary axis result 
with the Sommerfeld term subtracted. The result in Fig. 
5 is somewhat disappointing. Clearly, a linear depen- 
dence on temperature is obtained, and is easily traced to 
the linear dependence of Tt2ooiT) w —ttXT T >> uje 
(note that S2(T) — see Eq. (|17|l — has the same high 
temperature behaviour). Indeed, as x <C 1 in Eq. H22|l . 
one finds Koa{T) « K{0) - (2/7r)E2oo(r). However, the 
slope is evidently incorrect, and one may wonder how the 
discrepancy arises. 

A possibility presents itself when one examines more 
closely the imaginary part of the self-energy obtained an- 
alytically from Eq. ((njl and plotted in Fig. 6. The full 
expression at positive frequency is: 



Mz) = - 



TrXujf 



[coth(/3cjB/2) + f{z + uje) - f{z - lje)] , 



z < w/2 -UJE, 

Mz) = ^^[N{-LOE) + f{z-UJE)], 

W/2 -loe < z < W/2 + ue, 

S2(z) = 0, Z>W/2 + LUE. 



(23) 



(24) 
(25) 



As one can see, except for the small energy range below 
z = oje, ^2 given by Eq. (|23|l can be safely approximated 
by the Eq. H20|l above. However, near the band edge 
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Ko(T)(Eq. 13) 



KoJT) (Eq. 22) 



W = 1 eV 
Wp = 2 meV 



Kob(T) (Eq. 26) 
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z = W/2 —uje the self energy has a first step-hke change 
to a smaller value, and it then vanishes at z > W/2 + 
LUE- Viewed on the scale of Fig. 6, the fine structure 
present below oje and near the band edges is perhaps 
not important, while accounting for the upper cutoff z « 
W/2 clearly should be. This model can also be solved for 
analytically. The result is 



W 



f , 



[1 + a;^]-tan^^(-) 



1.2, 
— tan^V- 

TT X 



^ln(i- 



(26) 



where, as before, 



|S2oo(T)| 

W/2 



This result is shown in 



Fig. 5 as the blue dotted curve. Our attempt at improve- 
ment has resulted in a significantly inferior approxima- 
tion! 



FIG. 5: Various approximations of the kinetic energy (with- 
out the Sommerfeld expansion terms) compared with the nu- 
merical result (solid red curve). The green (dashed) curve is 
given by Eq. 122L where only the imaginary part of the self 
energy is included. It is not so far off the numerical result. 
However, the (relatively small) degree to which a discrepancy 
exists is serendipitous, as an obvious improvement (introduc- 
tion of a cutoff - see Fig. 6) results in the blue (dotted) curve. 
This latter 'improved' result is in far worse agreement with 
the numerics. 



T. 10 K 
T.50K 
T.100 K 
T=200 K 




-600 -400 -200 200 400 600 

w (meV) 

FIG. 6: Imaginary part of the self energy for various tempera- 
tures. At high temperatures, the use of a constant to describe 
the imaginary part of the self energy works very well; there is 
structure at the cutoffs barely visible on the scale shown here, 
but they are not very important. Despite this accurate ana- 
lytical representation of the imaginary part of the self energy, 
this approximation does not work well for the kinetic energy, 
as explained in Fig. 5. 



B. Real part of the self energy 

One can attempt to improve approximations to the 
imaginary part of the self energy (for example, by using 
a proper double step at the band edge), but this does 
not lead to any improvement in the subsequent Ko{T). 
Success is only attained when the real part of the self 
energy is included. A clear way to see this is as follows. 
One first notes that for an approximation which invokes 
a constant (with frequency) for the imaginary part of the 
self energy up to the bandwidth (i.e. S2(w -I- iS) = S200 
for < W/2, and zero otherwise), then Kramers-Kronig 
immediately implies that the real part of the self energy 
ReS = Si is given by a simple logarithm, i.e. 



lji[iLj + Id) = In 



W/2-UJ 



W/2 + UJ 



(27) 



Then the spectral function can be obtained analytically, 
and inserted into Eq. @, with the spectral function 
given in general by: 



I 



\^2{iO,T)\ 



T:{e~uj + Si (cj, T)f + Y.l{uj, T) ' 

The integral over e can be performed analytically, but the 
remaining integral over the frequency, w, must be done 
numerically. When this is done, the result for the kinetic 
energy is essentially indistinguishable from the numerical 
result at high temperatures. 

In light of the 'failures' of the approximations in the 
previous section, this requires further clarification. The 
(cij) frequency integral in Eq. extends from —00 to 
0. In the range —W/2 to zero the self energy has both 
a real and imaginary part; then the spectral function is 
some broad Lorentzian-like function. However, in the 
range —oo<uj< —W/2 the imaginary part of the self 
energy is zero, so the spectral function is given by a delta 
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FIG. 7: Real part of the self energy for various tempera- 
tures. Note that even at the lowest temperature shown (green, 
dashed curve), the overall slope at the origin looks positive 
on this scale; in fact there is some structure there, and the 
slope is in fact negative, but this plays almost no role in our 
results. This is because the most important contributions to 
the kinetic energy come from near the bottom of the band. 
Note that beyond the band edges the approximation given by 
Eq. (1281 is indistinguishable from the numerical result on this 
energy scale (a closer look would reveal some discrepancies — 
the numerical result has two logarithmic divergences at each 
band edge at low temperature, not one). 



function. The integral over the variable e is trivial, and 
the remaining integral over uj is non-zero, because the real 
part in this frequency range is also non-zero. If the real 
part of the self energy were zero, then this part of the 
integral (involving the delta function) is also zero, which 
is why we didn't mention it earlier with respect to the 
approximation Kqi,{T) (see Eq. (|26|) ). Here it must be 
included, and contributes a significant part (close to 20% 
of the total). 

This discussion firmly establishes the need to include 
the real part of the self energy in an accurate approx- 
imation for the kinetic energy. However, the previous 
exercise adds little to our fully numerical results, since 
one of the integrals must be done numerically. To see 
how well a simplified model does, we adopt the following 
simpler model. First, we take the imaginary part as be- 
fore, cut off at ±W/2. For the real part, if one expands 
Eq. (|2Z|) for \uj\ « W/2, the result is 



I]i(aj + iS) 



>.oo(T)|^ for 

|S2oc(r)l W/2 + LU 

In 



W/2 



\uj\ < W/2 
for 



w| > W/2. 
(28) 



In the first part, we have imposed a cutoff at the band 
edge. Note that the coefficient of uj is positive. This 



always occurs at high temperatures, consistent with the 
negative mass enhancement displayed in Fig. 4a. In fact, 
at high temperatures, |S2oo(r)| « nXT, and Eq. re- 
duces to the second term in Eq. H18|l . which dominates at 
high temperature, as Fig. 4 shows. Notice however that 
while approximation H28(l holds as soon a.s T ^ oje /2 the 
change in sign in Xo{T) occurs at higher temperature. 
In the second part, we have simply adopted Eq. 127|l . 
Fig. 7 shows how accurate this approximation is. In the 
frequency range \uj\ < W/2, the approximation is clearly 
crude; beyond this frequency range the approximation 
given by the second half of Eq. (|28|l is indistinguish- 
able from the numerical result. (Actually, if we use Eq. 
(|27|l . i.e. the logarithm for the entire frequency range, 
then this would be indistinguishable from the numeri- 
cal result over the entire frequency range. This accounts 
for our remarks above that insertion of Eq. H27() along 
with the constant imaginary part would result in a kinetic 
energy that is essentially indistinguishable from the full 
numerical result.) 

Explicitly, Eq. Q becomes 



2 

W 



— I duj 

w/2 



W/2 I-+W/2 

duj / de e6[uj — e — Si(cl; -|- iS)] + 

o J-W/2 
+W/2 ^ 

de e— - — 

(e ■ 



j2oo 



{T)\ 



^-HSi(c. + z<5))2-HS2^(r)- 



(29) 



In the first line, the e integration is trivial. Care must be 
taken concerning the bandedge cutoffs, so that the uj inte- 
gration no longer extends to — cx3. This lower cutoff, Wc, is 



determined by the condition — 



In 



W/2+uj^ 
W/2-i^a 



—W/2, which can easily be solved by a few iterations. 
The remaining lu integral is elementary; we write it as 



Kl{T) 



W 



2 

1 + -X 
TT 



In 



V 



I' 



yln 



(30) 

where x is as before, and y = ^j^. The ratio y is always 
greater than unity. 

In the second line we write Si — —Xlu, where X — ——x 
represents a negative mass enhancement. Substitution 
into the integral allows this piece (which we denote 
i^Q^(T)) to be evaluated analytically as before: 



kUt) 



-(1 + A)ln 

TT 

(A(A-f 2) ~x 



W 1 

4 l + X 



(2 + A)2 + 



(l + a;2)2tan-i(l) 



.2^iftan-i(^)-tan-i(- 

TT V X X' 



(31) 

with X = ' " . The approximation Koc{T) consists of 
the sum of these two pieces, Koc{T) = K^ci^) + Kl^{T). 
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FIG. 8: Comparison of Koc{T) approximation (solid green 
circles) with numerical results (solid red curve). Agreement 
is excellent. Also shown are the results for Kq^{T) alone, i.e. 
without accounting for the delta-function part. Clearly both 
real and imaginary parts of the self energy are necessary for 
accurate results. 



FIG. 9: Comparison of analytical model [KodT)) (symbols) 
with numerical results (curves) for several boson frequencies. 
Not the expanded vertical scale. Clearly, the model deteri- 
orates as the boson frequency increases (also note that the 
result is also no longer linear). Nonetheless, the discrepancy 
for uJ_B = 25 meV is still less than 2%. 



IV. DISCUSSION AND SUMMARY 



These are plotted for cj^ = 2 mcV and A = 1 in Fig. 
8 (solid green circles). Agreement with the numerical 
result (solid red curve) is now excellent, especially when 
one considers the crudeness of the approximation for the 
real part of the self energy shown in Fig. 7. We also 
show the result for Kq^{T), i.e. without the correction 
from the delta-function piece (which arises only because 
the real part of the self energy is non-zero) with the open 
blue circles. Clearly both real and imaginary parts of the 
self energy are required for accurate results. 

Finally, we show results in Fig. 9 for several boson 
frequencies, over a somewhat higher temperature range. 
The approximations used to obtain Kqc{T) clearly dete- 
riorate as the boson frequency increases with respect to 
the temperature scale involved. Note in particular that 
in the temperature regime where the result is no longer 
linear in temperature the qualitative agreement is not so 
good (i.e. off by more than a percent). It is also ap- 
parent in this figure that the crude approximation for 
the real part of the self energy (first part of Eq. (|28|) ) 
results in some minor discrepancies at higher tempera- 
tures. As mentioned before, replacement of this linear 
approximation with a single logarithm (still an approxi- 
mate result) makes the 'analytic' results agree essentially 
perfectly with the numerical results. 



In the previous sections we have described the effects of 
the interaction of electrons with a soft collective mode on 
the electronic kinetic energy, which is in turn related to 
the optical sum rule probed by the experiments. In con- 
trast to the general expectation for Fermi- liquid systems, 
here the main contribution to the temperature variation 
of the kinetic energy cannot be attributed to the Som- 
merfeld term, i.e. to the thermal smearing of the Fermi 
occupation number. Instead, the relevant effects come 
from the temperature dependence set by the interaction 
itself, via the electronic self energy. While discussing 
the sum rule behavior the presence of a finite bandwidth 
W is a, crucial ingredient, since the temperature varia- 
tions are expected to vary as a power of T/W or £2/!^. 
The latter case is relevant here, with the kinetic energy, 
as given by Eq.s (|30|) - (|31|l . scaling as Yj2/W . Moreover, 
since Yi2{T) — —nXT at high temperatures, we also found 
a linear temperature depletion of the optical sum. How- 
ever, in deriving this result we had to properly take into 
account also the structure of the real part of the self en- 
ergy, whose behavior at low frequency is quite strongly 
affected by the presence of a finite-bandwidth cut-off, in 
contrast to the general expectation coming from the anal- 
ysis of the infinite-bandwidth approximationi^iiSS 

Finally, we would like to present a simple argument 
which can help the reader to understand the difference 
between the effect of thermal smearing of the Fermi func- 
tion with respect to the effect of the interaction consid- 
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ered here. Going back to the expression (0) of the kinetic 
energy, we can understand the Sommerfeld term as due to 
the fact that by varying the temperature one has a vari- 
ation of the occupation factor Sn{(_) of order f/T over a 
range ~ T around the Fermi level (here ^ = e — fi) and 
otherwise so that: 



SK 

W 



I j-T T"^ ( T 



(32) 



As a consequence, the power is due to the fact that 
the integration of energy, which gives e^, on a range 
T jW ^ produces a power of (T jWY in the kinetic en- 
ergy. The question then arises: do interactions simi- 
larly simply change the occupation number in a region 
I]2|/M^ around the Fermi level, which would then lead 
to a (|S]2|/VF)^ dependence of the kinetic energy for the 
interacting case? The answer we have found is clearly, 
no. Indeed, when the occupation number coming from 
Eq. © is considered, using the same approximations 1)20(1 
and (|28l) considered above for the self-energy we can esti- 
mate no(e) (i.e. the contribution coming from the T = 
limit of the Fermi function, in analogy with the definition 
® of isTolr)) as: 



"-o(e) = 



1 



^(1 + A) 



tan 



.1 t^(l + A)/2 + , 



— tan 



■J2oo| 



IS200I 

(33)- 



As a consequence, one can roughly say that no(e) ~ 1 — 
Q!|S2oo|/VF, with some coefhcient a of order 1, in most of 
the energy range between — Wjl and 0, and not only in a 
small region around the Fermi level. In other words, the 
interaction affects the occupation number even far from 



the Fermi level, so that the energy integration in Eq. @ 
is now scaling as: 



W 



W/2 



W/2 



aW 



W 



(34) 



Thus, the kinetic energy is directly proportional to the 
temperature variation of the imaginary part of the self 
energy, even though the exact coefficient a should be de- 
termined taking care also of contributions from the real 
part of the self-energy. Even though the previous rela- 
tion is not valid for stronger interactions (as in Hubbard- 
like modelsi^) , one can reasonably expect that by model- 
ing the inelastic scattering with a more general phonon- 
like spectrum, the sum rule behavior will roughly fol- 
low the modification introduced in the electronic self- 
energy in the way indicated by Eq. . This possibility 
leads to an understanding of the experimental observa- 
tion of a non-quadratic temperature dependence of the 
sum rule in some cuprate superconductors, where also 
the high-energy scattering rate shows a linear tempera- 
ture scaling^. 
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